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SUMMARY 


The  work  summarized  in  this  report  covers  the  seventh  quarter  of  a  program 
with  a  goal  that  is  threefold:  first,  to  study  the  properties  of  native  point  defects  in 
infrared  focal-plane  array  (IRFPA)  active  and  substrate  materials;  second,  to  study 
the  properties  of  native  point  defects  in  two  classes  of  photonic  materials,  the  wide-gap 
II- VI  compounds  (ZnSe  as  the  prototype  for  which  impurity  properties  will  also  be 
calculated)  and  the  nonlinear  optical  materials  (LiNbOs  as  the  prototype);  and  third 
to  study  the  properties  of  HgZnTe  as  a  very-long-wave  infrared  (VLWIR)  detector 
material.  Our  accomplishments  in  the  seventh  quarter  include 

•  Beginning  of  calculation  of  the  tellurium  antisite  -  mercury  vacancy 
pair.  Prehminary  results  indicate  that  these  complex  concentrations 
will  be  high  enough  to  impact  device  performance.  They  serve  as 
reservoirs  of  Tejjg ,  which  keeps  the  concentrations  of  these  donors  high 
in  commonly  used  annealing  procedures,  and  the  complexes  plus  the 
TeHg  may  account  for  the  Shockley-Read  centers  that  currently  limit 
FPA  performance  at  intermediate  operating  temperatures.  Annealing 
procedures  to  eliminate  the  adverse  effects  of  these  defects  are  being 
formulated. 

•  Continued  calculation  of  the  defect  ionization  energies  in  CdTe,  ZnSe, 
and  LiNbOa. 

•  Calculation  of  the  temperature  dependent  band  gap  of  semiconductors. 

Our  calculation  predicts  the  temperature  variations  of  the  band  gaps 
and  band  offsets  of  HgCdTe,  GaAs,  InSb,  InAs,  and  InP  correctly.  The 
band  offset  between  HgTe  and  CdTe  varies  from  350  meV  at  T=0  to 
215  meV  at  300  K.  This  will  modify  the  heteroj unction  designs  for 
devices  intended  to  operate  at  elevated  temperatures. 
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1.  NATIVE  POINT  DEFECTS  IN  HgCdTe  AND  RELATED  IR  MATERIALS 

1.1  MERCURY  VACANCY  TELLURIUM  ANTISITE  COMPLEX 

We  have  continued  to  calculate  the  formation  energy  for  the  mercury  vacancy 
tellurium  antisite  pair.  The  revised  calculations  will  include  the  relaxation  of  the  two 
defects.  Once  the  formation  energy  is  calculated,  we  will  revise  our  estimates  of  the 
density  of  this  defect  and  its  impact  on  device  properties. 

1.2  NATIVE  POINT  DEFECTS  IN  CdTe 

We  worked  on  a  method  for  calculating  the  ionization  energies  of  the  defects  in 
semiconductors.  While  we  were  able  to  assume  that  the  ionization  energies  of  donors 
and  acceptors  are  zero  in  HgCdTe  because  of  its  small  band  gap,  we  cannot  do  this 
for  CdTe. 

1.3  TEMPERATURE-DEPENDENT  BAND  GAP 

The  electron-phonon  interaction-induced  band-edge  shifts  in  semiconductors  axe 
calculated  using  accurate  band  structures.  The  temperature  variation  of  gaps  in  GaAs 
(done  to  validate  our  method)  and  Hgo.78Cdo.22Te  (MCT)  have  been  calculated  and 
are  found  to  agree  well  with  experiments.  While  the  simple  picture  that  the  intra 
(inter)-band  transitions  reduce  (increase)  the  gap  still  holds,  we  have  shown  that 
both  the  conduction  band  edge  Ec  and  valence  band  edge  move  down  in  energy. 
These  large  shifts  in  E^  can  affect  the  valence  band  offsets  in  hetero junctions  at  finite 
temperature.  Detailed  band-to-band  analysis  explains  why  the  gap  in  InSb  decreases 
with  T,  whereas  it  increases  in  MCT  alloys  with  a  comparable  band  gap.  A  paper 
on  this  work  (included  here  as  the  Appendix)  has  been  submitted  to  Physical  Review 
Letters. 

2.  HgZnTe  AS  AN  LWIR  DETECTOR  MATERIAL 

We  have  completed  the  calculation  of  the  band  structure  of  Hgi_iZni,Te  as  a 
function  of  x  using  a  hybrid  pseudopotential  tight-binding  method.  This  method  has 
been  applied  extensively  to  a  wide  range  of  group  IV,  IITV,  and  IT VI  compounds 
and  alloys!  with  a  great  degree  of  success.  The  method  uses  a  generic  pseudopotential 
for  all  of  the  semiconductors  to  which  a  tight-binding  Hamiltonian  is  added  as  a  per¬ 
turbation.  The  parameters  of  the  tight-binding  Hamiltonian  are  chosen  empirically 
so  as  to  give  good  agreement  with  the  experimental  band  structure.  The  coherent 
potential  approximation  (CPA)  is  used  to  calculate  the  alloy  band  structure  from 
those  of  the  constituent  compounds. 


!  A.-B.  Chen  and  A.  Sher,  “The  theory  of  semiconductor  alloys”,  unpublished. 
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3.  WIDE-GAP  ll-VI  COMPOUNDS  (ZnSe  AS  THE  PROTOTYPE) 

Our  work  has  continued  on  the  calculation  of  the  ionization  energies  of  the  defects 
in  ZnSe. 

4.  NONLINEAR  OPTICAL  MATERIALS  (LINbOs  AS  THE  PROTOTYPE) 

This  quarter  our  work  has  continued  on  the  calculation  of  the  ionization  energies 
of  the  defects  in  LiNbOa. 

5.  WORK  PLANNED 

In  the  next  quarter  we  will  continue  the  calculation  of  the  ionization  energies  of 
the  defects  in  CdTe,  ZnSe,  and  LiNbOa,  and  the  incorporation  of  these  energies  into 
the  quasichemical  formalism  to  predict  densities  of  the  ionized  defect  concentrations. 
We  expect  to  complete  the  calculation  of  the  binding  energy  of  the  tellurium  antisite 
mercury  vacancy  pair.  We  will  be  calculating  the  transport  properties  of  HgZnTe, 
including  the  electron  and  hole  mobilities  and  the  Hall  coefficient  and  carrier  concen¬ 
tration  as  a  function  of  x,  T,  Nc,  and 
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APPENDIX 


Temperature  dependence  of  band  gaps  in  semiconductors 


Srinivasan  Krishnamurthy®,  A.-B.  Chen\  and  A.  Sher® 

“  SRI  International,  Menlo  Park,  CA  94025 
’’  Physics  Department,  Auburn  University,  Auburn,  AL  S6349 
(May  25,  1994) 


The  electron-phonon  interaction-induced  band-edge  shifts  in  semiconductors  are  calculated 
using  accurate  band  structures.  The  calculated  temperature  variation  of  gaps  in  GaAs  and 
IIgojsCdo.22Te{M.CT)  agrees  well  with  experiments.  While  the  simple  picture  that  the  intra  (inter)- 
band  transitions  reduce  (increase)  the  gap  still  holds,  we  show  that  both  the  conduction  band  edge 
Ec  and  valence  band  edge  Ev  move  down  in  energy.  These  large  shifts  in  Ev  can  affect  the  valence 
band  offsets  in  heterojunctions  at  finite  temperature.  Detailed  band-to-band  analysis  explains  why 
the  gap  in  InSb  decreases  with  T,  whereas  it  increases  in  MCT  alloys  with  a  comparable  band  gap. 


PACS  numbers  63.20.Kr,  71.25.Tn 

The  temperature  (T)  dependence  of  energy  gaps  of 
semiconductors  is  of  great  physical  and  technological 
interest.  Numerous  theoretical  [1-9]  and  experimental 
[10-19]  studies  have  been  undertaken  to  obtain  both 
qualitative  and  quantitative  variations  of  various  gaps 
in  semiconductors.  While  the  gap  decreases  with  in¬ 
creasing  temperature  in  medium-gap  and  wide-gap  semi¬ 
conductors,  it  increases  in  small-gap  materials  such  m 
HgCdTe,  PbS,  PbSe,  and  PbTe.  The  temperature  vari¬ 
ation  is  usually  attributed  to  two  causes:  thermal  ex¬ 
pansion  of  the  lattices  and  electron-phonon  interactions. 
Thermal  expansion  always  reduces  gaps.  Here  we  show, 
starting  from  accurate  band  structures  and  wave  func¬ 
tions,  from  proper  phonon  dispersion  relations,  and  tak¬ 
ing  account  of  matrix  elements  of  the  electron-phonon 
interactions,  that  the  calculated  gap  variations  in  GaAs 
and  B'go.7aCda.22Te  agree  well  with  experiments.  In 
a  perturbation-theory  treatment  of  electron-phonon  in¬ 
teractions,  the  intrarband  transitions  reduce  the  gap 
whereas  interband  transitions  increase  it,  and  the  net 
shift  in  the  gap  can  be  positive  or  negative.  We  show 
that  while  this  simple  physical  picture  is  still  valid,  both 
conduction  and  valence  band  edges  move  down  in  energy. 
When  the  valence  band  edge  moves  more  than  the  con¬ 
duction  band  edge,  the  gap  increases  with  T  as  in  the 
case  of  Hgo.7aCdQ,22Te.  The  reverse  occurs  for  GaAs 
and  most  other  semiconductors.  This  observation  has  an 
important  effect  on  our  understanding  of  the  variation 
of  band  offsets  in  semiconductor  heterojunctions.  The 
contributions  &om  each  phonon  branch  to  each  electron 
band  have  been  obtained  to  assist  physical  understanding 
of  the  underlying  causes  of  the  variations. 

Our  calculation  of  the  temperature  dependence  of  the 
band  gap  starts  with  accurate  band  structures.  Empiricad 
pseudopotential  form  factors  are  used  to  construct  a  hy¬ 
brid  tight-binding  (HTB)  Hamiltonian  in  a  minimum  set 
of  sp^  Slater  orbitals  per  atom.  This  Heimiltonian  is  then 
transformed  into  an  orthonormal  basis.  A  site-diagonal 


spin-orbit  Hamiltonian  and  a  perturbative  Hamiltonian 
are  added  to  fine-tune  the  band  structures  to  agree  well 
with  experiments  [20,21].  Various  results  obtained  us¬ 
ing  these  band  structures  are  found  to  be  quite  reliable 
[21-23].  This  treatment  tests  the  accuracy  of  the  wave 
functions  as  well  as  the  energies. 

The  dilation  contribution  to  the  band  gap  reduction 
is  given  [5,9]  by  3  aj  B  dEg/d  P,  where  the  thermal 
expansion  coefficient  of  the  lattice  ot,  the  bulk  modu¬ 
lus  B,  and  the  change  in  the  gap  with  pressure  are  ob¬ 
tained  from  the  literature  [19].  The  electron-phonon  in¬ 
teractions  that  cause  the  band  structures  to  change  are 
treated  in  perturbation  theory.  The  total  potential  is  as¬ 
sumed  to  be  a  sum  of  potentials  from  single  atoms.  The 
atomic  potential  in  the  solid  is  traditionally  expanded 
in  a  Taylor  series  with  only  the  leading  term  retained, 
Sind  the  energy  is  evaluated  in  a  second-order  perturba¬ 
tion  theory  commonly  known  as  a  self  energy  correction. 
However,  it  has  been  demonstrated  by  a  number  of  re¬ 
searchers  [3,4,6]  that  retention  of  first-order  perturbation 
terms  with  a  second  term  in  Taylor  series  is  necessary  to 
preserve  symmetry.  We  retain  both  terms.  The  change 
in  the  energy  at  a  given  wave  vector  k  is 


AEnit  =<  nk|V2|nk  >  + 

n'k' 


<n'k'|Vi|nk>  p 

-  -^n'k' 


(1) 


where  Vi  and  V2  are  the  first  two  terms  in  the  Tay¬ 
lor  expansion  of  the  total  potential  in  powers  of  atomic 
displacements  In  the  TB  formalism,  Eq.  (1)  can  be 
written  in  terms  of 

<  Vj'a'\Vi\\ja  >=  VV„„-(df/;')  •  -  $y),  (2) 

and 

<  Vj'a'mja  >  =  •  (V)Vaa'(d{/;')  • 

+6;  ■  (V)^  Vaa'C*^!/'  ) 
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t  •  .  / 

where  dj//  is  the  position  vector  connecting  atomic  sites 
1,  species  (anion  or  cation)  j  and  site  1',  species  j',  and 
Vaa'  (di/;  )  is  a  HTB  matrix  element  between  the  orbitals 
a  and  a'  located  on  those  atoms.  From  quantum  theory 
of  harmonic  crystals,  the  atomic  displacements  ^  can  be 
expressed  in  terms  of  normal  modes  of  phonons.  We  have 

fu  =  [57m7l‘ 

(-1) 

where  q  and  u  are  phonon  wave  vector  and  frequency,  A 
denotes  phonon  branch,  a(at)  is  the  phonon  annihilation 
(creation)  operator,  M  is  the  atomic  mass  and  e  is  an 
eigenvector  of  the  six-dimensional  eigenvalue  problem; 

Mu^e  =  D(q)e  (5) 

Evaluation  of  the  matrix  elements  given  by  Eqs.  (2) 
and  (3)  requires  knowledge  of  spatial  variations  of  the 
interatomic  TB  matrix  elements  Vaa'-  In  Harrison’s  uni¬ 
versal  TB  approach  [24],  these  matrix  elements  scale  as 
d“^.  In  our  generalization,  we  assume  that  Vaa'  varies  as 
d“'"  and  the  repulsive  first-neighbor  pair  energy,  follow¬ 
ing  Harrison’s  overlap  argument,  as  Tj/d^"*.  The  two  un¬ 
knowns  m  and  t)  are  determined  by  requiring  that  the  cal¬ 
culated  equilibrium  bond  length  and  bulk  modulus  agree 
well  with  experiments.  This  approach,  with  electrons 
and  phonons  treated  from  the  same  underlying  Hamil¬ 
tonian,  hcis  previously  been  used  successfully  to  explain 
hot  electron  transistor  characteristics  [23]  and  is  also  in 
fairly  good  agreement  with  first-principles  calculations 
[25].  The  dynamical  matrix  D  is  calculated  from  the 
valence  force  field  model  [26]. 

The  calculational  procedure  is  as  follows.  For  a  chosen 
material,  m  and  t]  are  evaluated.  Then  the  first  and 
second  derivatives  of  all  interatomic  matrix  elements  are 


FIG.  1.  Change  in  the  band  gap  of  GaAs  with  temperature. 


obtained.  The  dynamical  matrix  is  diagonalized  to  ob¬ 
tain  u  and  e  as  a  function  of  q  and  A.  The  phonon 
structures  and  electronic  band  structures  are  used  [Eqs. 
(1)  through  (4)]  to  obtain  the  change  in  the  band  energy 
at  a  given  k.  The  polar  coupling  terms  are  included  in 
the  longitudinal  optical  phonon  contributions.  As  we  are 
interested  in  studying  the  change  in  the  direct  gap,  k  is 
taken  to  be  zero. 

The  calculated  band-gap  change  as  a  function  of  T  in 
GaAs  and  Hgo.7sCdo.22Te  is  shown  in  Figs.  1  and  2, 
respectively.  With  increasing  T,  the  direct  gap  decreases 
in  GaAs  but  increases  in  Hgo.7sCdo.22Te.  Notice  that, 
in  both  cases,  the  calculated  values  aje  typically  within 
10  to  15  meV  of  experimental  values  [10,13,16-18].  The 
cross  (x)  at  T=0  represents  the  calculated  zero-point 
correction  to  the  gap  (-26  meV  for  GaAs  and  10  meV 
for  Hgo.7sCdo.22Tt).  The  zero  temperature  band  gap 
calculated  without  electron-phonon  interactions  should 
have  this  correction  subtracted  for  comparison  them  to 
experimental  values. 


FIG.  2.  Change  in  the  band  gap  of  Hgo.TtCdo.33Te  with 
temperature. 

The  change  in  the  gap  is  traditionally  explained  in 
terms  of  inter-  and  intraband  interactions.  The  intrava¬ 
lence  (conduction)  bemd  interactions  push  the  valence 
(conduction)  band  edge  up  (down),  thus  reducing  the 
gap.  Similarly,  the  valence-conduction  band  interactions 
increase  the  gap.  Hence  one  might  expect  the  gap  to  de¬ 
crease  in  wide-gap  semiconductors  and  possibly  increase 
in  small-gap  semiconductors.  In  addition,  arguments 
based  only  on  total  density  of  states  and  ignoring  varia¬ 
tions  in  matrix  elements  will  predict  both  the  conduction 
band  edge  £?«  and  the  valence  band  edge  move  up  in 
energy,  particularly  in  small-gap  materials,  because  the 
hole  effective  mass  is  one  to  two  orders  of  magnitude 
larger  than  the  electron  mass.  As  seen  from  Fig.  3,  our 


A-2 


‘fig.  3.  Variation  of  conduction  and  valence  band  edges 
(solid  lines)  with  temperature.  Interaction  with  valence  bands 
pushes  the  states  up  (dashed  lines)  and  that  with  conduction 
bands  pushes  them  down  (dotted  lines). 


detailed  calculations  do  not  support  this  traditional  view. 
We  find  that  both  Eg  and  (solid  lines)  move  down 
in  energy.  This  trend  is  observed  in  all  semiconductors 
tried  (InP,  InAs,  InSb,  HgCdTe  alloys,  InTlP  alloys,  and 
GaSb).  The  movement  of  the  band  edges  due  to  inter¬ 
action  with  valence  bands  (dashed  lines)  and  that  with 
conduction  bands  (dotted  lines)  is  also  shown  in  Fig.  3. 
We  see  that  interaction  of  the  band  edges  with  the  con¬ 
duction  bands  (CBs)  is  much  stronger  than  with  valence 
bands  (VBs),  and  consequently  both  E^  and  Ec  move 
down  in  energy. 

TABLE  I.  Calculated  change  at  300  K  in  the  valence  (v) 
and  conduction  (c)  band  edge  energies  in  meV.  Contributions 
from  interaction  with  various  phonon  modes  are  shown  in  the 
third  through  tenth  rows. 


band 

1 

2 

3 

4 

5 

6 

7 

8 

V 

0.3 

3.0 

9.4 

9.7 

-36.4 

-37.9 

-42.5 

-30.3 

Total 

c 

0.7 

1.8 

13.1 

18.6 

-97.0 

-33.8 

-68.3 

-36.2 

V 

0.2 

1.8 

6.8 

4.8 

-23.7 

-27.3 

-30.0 

-15.9 

TA 

c 

0.4 

0.6 

11.0 

16.4 

-11.5 

-15.0 

-52.0 

-14.2 

V 

0.0 

0.6 

0.4 

1.5 

-6.9 

-3.1 

-2.7 

-5.4 

LA 

c 

0.2 

0.7 

0.7 

0.6 

-62.2 

-10.1 

-1.7 

-8.5 

V 

0.0 

0.4 

1.8 

2.9 

-3.3 

-6.3 

-5.5 

-3.8 

TO 

c 

0.0 

0.1 

1.1 

1.4 

-5.4 

-5.2 

-12.9 

-9.2 

V 

0.0 

0.2 

0.5 

0.6 

-2.5 

-1.2 

-4.4 

-5.3 

LO 

c 

0.1 

0.3 

0.3 

0.2 

-18.0 

-3.5 

-1.7 

-4.3 

We  analyze  the  strength  of  inter-  and  intraband 
electron-phonon  interactions  by  presenting  the  contribu¬ 
tions  from  each  band  and  from  each  phonon  mode.  Table 
I  lists  the  calculated  values  for  GaAs  at  300  K.  Although 
spin  is  included  in  our  band  structure  calculations,  only 
spin  averaged  values  are  listed.  Band  indices  1  to  4  cor¬ 
respond  to  VBs,  and  the  others  to  CBs.  The  first  two 
rows  show  changes  in  E^  and  Ee  due  to  interactions  with 
various  bands.  Contributions  firom  each  phonon  mode 
are  listed  in  the  remaining  rows.  The  lowest  VB  is  about 
13  to  15  eV  below  E^  and  Ec,  and  hence  the  interac¬ 
tion  impacts  the  beind  edges  only  very  little.  This  band 
could  be  omitted  without  much  loss  of  accuracy.  We 
see  that  the  interaction  with  other  VBs  is  also  generally 
weak  and  the  matrix  elements  connecting  VB  states  are 
always  found  to  be  small.  However,  interaction  of  the 
beind  edges  with  CBs  is  very  strong  emd  negative.  Pm- 
ticularly,  the  interaction  of  Eg  with  the  lowest  conduction 
band  is  strong  because  of  a  combination  of  large  matrix 
elements  and  small  energy  denominators  [Eq.  (1)].  As  Eg 
moves  down  more  than  E^ ,  the  gap  narrows  in  GaAs.  To 
understand  the  role  of  various  phonons,  we  display  the 
contribution  from  each  mode  separately.  The  calculated 
values  of  change  in  the  band  edges  at  300  K  in  GaAs  are 
listed  in  the  third  through  tenth  rows  of  the  table.  While 
the  conduction  band  edge-CB  interaction  mediated  by 
optical  phonons  contribute  modestly  to  the  value  of  Eg, 
we  see  that  acoustic  phonons  account  for  about  75%  of 
the  total  change  in  the  valence  and  conduction  band  edge 
energies. 

Our  calculations  for  Hgo.78Cdo,22Te  alloy  show  a  qual¬ 
itatively  similar  role  for  phonons.  However,  we  find  that 
Ey  moves  down  more,  because  of  the  narrow-gap  (100 
meV),  than  Eg  and  that  the  gap  increases  with  temper¬ 
ature.  Then  we  treated  InSb  and  Hgo.70Cdo.30Te.  Al¬ 
though  they  both  have  the  same  zero  temperature  gap 
of  0.235  eV,  the  InSb  gap  decreases  with  T,  whereas 
the  Hgo.7oCdo.3oTe  gap  increases  with  T.  The  difference 
is  mainly  because  of  the  lowest  conduction  band  width 
and  the  intervalley  separation.  In  the  HCT  eJloy,  the  F 
and  the  next  nearest  valley  L  are  separated  by  2.7  eV 
whereas  that  energy  separation  is  only  0.7  eV  in  InSb. 
The  smaller  energy  denominator  [Eq.  (1)]  in  the  inter¬ 
action  of  the  conduction  band  edge  with  the  lowest  con¬ 
duction  band  pushes  down  Eg  much  more  in  InSb  than 
in  Hgo.7oCdo.3oTe,  and  consequently  the  gap  decreases 
in  InSb.  Of  course,  for  a  larger  Cd  concentration  the 
zero  T  gap  is  also  larger,  and  E„  moves  down  by  only 
a  small  amount  and  the  finite  temperature  gap  starts  to 
decrease.  The  combination  of  gap  size,  the  lowest  con¬ 
duction  band  width,  and  intervalley  separation  gives  rise 
to  these  interesting  variations  in  the  gap  with  T. 

The  observation  that  both  Ey  and  Eg  move  down 
in  energy  has  an  important  effect  on  band  offsets  in 
heterojunction-based  devices.  For  example,  the  zero 
temperature  valence  band  offset  between  iffl'o.zsC'do.zzT^e 
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and  CdTe  is  believed  to  be  around  350  meV.  However,  we 
find  that  at  300  K,  in  HgQ,7sCdo,22l'^  S'Hd  in  CdTe 
moves  down  by  215  meV  and  30  meV,  respectively.  If 
the  dipole  contribution  remains  the  same,  the  valence 
band  offset  decreases  to  165  meV  at  300  K.  In  addition 
to  the  electron-phonon  interactions  discussed  above,  lat¬ 
tice  dilation  changes  the  band  edges  differently  [9].  In 
any  case,  this  band  offset  change  has  enormous  influence 
on  the  design  of  IR  absorption  and  l^ing  devices. 

In  addition  to  GaAs  and  ifiro.78CdQ.22Te,  we  studied 
the  band-gap  variation  with  temperature  in  InAs,  InP, 
InSb,  GaSb,  and  CdTe  compounds.  The  calculations 
produced  correct  trends  in  all  these  compounds,  and  the 
values  in  InP,  GaSb,  and  CdTe  are  in  fairly  good  agree¬ 
ment  with  experiments.  However,  the  calculated  changes 
in  the  bamd  gap  of  InAs  and  InSb  were  about  a  factor  of 
2  smaller  than  in  the  experiments.  We  And  that  our 
calculated  TA  phonon  frequencies  away  from  zone  cen¬ 
ter  in  these  compounds  were  considerably  larger  than  in 
experiments.  As  noted  from  Table  I,  a  substantial  contri¬ 
bution  comes  from  acoustic  phonons.  Consequently,  our 
theoretical  values  of  Ej(T)  are  smaller  than  in  experi¬ 
ments.  Although  these  inaccuracies  do  not  change  the 
conclusions  of  this  Letter,  better  predictability  should 
result  from  improvement  in  the  dynamical  matrix  cal¬ 
culated  from  the  underlying  Hamiltonian.  In  addition, 
at  higher  temperatures  higher-order  perturbation  terms 
must  be  included  along  with  finite  temperature  ‘renor- 
maUzed’  bands  rather  than  the  zero  temperature  bands. 
Such  renormalization  affects  the  monotonic  change  in  the 
gap  and  introduces  nonlinear  terms. 

In  summary,  we  have  calculated  temperature  varia¬ 
tions  of  band  gaps  in  veirious  semiconductors.  A  fairly 
accurate  HTB  Hamiltonian  is  used  in  the  calculation  of 
electron  and  phonon  structures.  The  calculations  explain 
the  increase  in  the  band  gap  of  iffl'o.zsC'do.zzTe,  and  the 
decrease  in  the  band  gap  of  GaAs,  and  show  that  acous¬ 
tic  phonons  make  the  major  contribution.  Contreiry  to 
traditional  thinking  based  on  total  density  of  states  ar¬ 
guments,  we  find  that  both  valence  and  conduction  band 
edges  move  down  in  energy.  One  important  consequence 
of  this  observation  will  be  in  the  band  offsets  in  semicon¬ 
ductor  heterojunction  devices.  Finally,  there  is  a  small 
and  always  negligible  zero-point  motion  contribution  to 
low-temperature  band  gaps  arising  from  electron-phonon 
interactions. 

We  thank  Dr.  M.  Cardona  of  the  Max  Planck  In¬ 
stitute,  Stuttgart,  for  pointing  us  to  several  references. 
Funding  from  ONR  (contract  N00014-93-C-0091)  and 
ARPA  (contract  MDA972-92-C-0053)  is  gratefully  ac¬ 
knowledged. 


[1]  Y.  P.  Varshini,  Physica  34,  149  (1967). 

[2]  V.  Heine  and  J  .A.  Van  Vechten,  Phys.  Rev.  B13,  1622 
(1976). 

[3]  P.  B.  AUen  and  V.  Heine,  J.  Phys.  C  9,  2305  (1976) 

[4]  P.  B.  Allen  and  M.  Cardona,  Phys.  Rev.  B  27,  4760  (1983) 

[5]  S.  Gopalan,  P.  Lautenschlager,  and  M.  Cardona,  Phys. 
Rev.  B  36,  5577  (1987) 

[6]  M.  Cardona  and  S.  Gopalan,  in  Progress  on  electron  prop¬ 
erties  of  solids,  R.  Girlanda  et  al.,  eds.  (Kluwer  Academic 
publishers,  1989),  p.  52. 

[7]  R.  D.  King-Smith,  R.  J.  Needs,  V.  Heine,  and  M.  J.  Hodg¬ 
son,  Europhys.  Lett.  10,  569  (1989) 

[8]  S.  ZoUner,  S.  Gopalan,  and  M.  Cardona,  Sol.  State  Comm. 
77,  485  (1991) 

[9]  K.  J.  Malloy  and  J.  A.  Van  Vechten,  J.  Vac.  Sci.  Technol. 
B  9,  2112  (1991). 

[10]  P.  Lautenschlager,  M.  Garriga,  S.  Logothetidis,  and 
M.  Cardona,  Phys.  Rev.  B  35,  9174  (1987);  and  references 
cited  therein. 

[11]  Z.  Hang,  D.  Yan,  F.  H.  Poliak,  G.  D.  Pettit,  and 
M.  Woodall,  Phys.  Rev.  B  44,  10546  (1991). 

[12]  L.  Pavisi,  F.  Piazza,  A.  Rudra,  J.  F.  Carlin,  and 
M.  Ilegems,  Phys.  Rev.  B  44,  9052  (1991). 

[13]  E.  Grilli,  M.  Guzzi,  R.  Zamboni,  and  L.  Pavesi,  Phys. 
Rev.  B  45,  1638  (1992). 

[14]  P.  Y.  Liu  and  J.  C.  Maan,  Phys.  Rev.  B  47, 16274  (1993). 

[15]  M.  E.  AUali,  C.  B.  Sorenson,  E.  Veje,  and  P.  T.  Petersson, 
Phys.  Rev.  B  48,  4398  (1993). 

[16]  G.  L.  Hansen,  J.  L.  Schmit,  and  T.  N.  Casselman,  J. 
Appl.  Phys.  63,  7099  (1982). 

[17]  D.  G.  Seiler,  J.  R.Lowney,  C.  L.  Littler,  and  M.  R.  Loloee, 
J.  Vac.  Sci.  Technol.  A  8,  1237  (1990). 

[18]  J.  C.  Brice,  Properties  of  mercury  cadmium  telluride, 
J.  Brice  and  P.  Capper,  eds.  (EMIS  datareviews  series  3, 
INSPEC  publication,  1987),  p.  105. 

[19]  Landolt-Bornstein  Numerical  data  and  functional  rela¬ 
tionships  in  science  and  technology,  Madelung,  Schultz  and 
Weiss  (eds.),  New  series,  Vol.  17  (1982). 

[20]  A.  -B.  Chen  and  A.  Sher,  Phys.  Rev.  B  23,  5360  (1981). 

[21]  S.  Krishnamurthy,  A.  -B.  Chen,  and  A.  Sher,  J.  Appl. 
Phys.  63,  4540  (1988). 

[22]  S.  Krishnamurthy,  A.  Sher,  and  A.  -B.  Chen,  Phys.  Rev. 
Lett.,  64,  2531  (1990);  Appl.  Phys.  Lett  65,  1002  (1989); 
Appl.  Phys.  Lett  62,  468  (1988). 

[23]  S.  Krishnamurthy,  A.  Sher,  and  A.  -B.  Chen,  Appl.  Phys. 
Lett.  53,  1853  (1988). 

[24]  W.  Harrison,  Electronic  structure  and  properties  of  solids 
(Freeman,  San  Francisco,  1980). 

[25]  S.  Krishnamurthy  and  M.  Cardona,  J.  Appl.  Phys.  74, 
2117  (1993). 

[26]  R.  M.  Martin,  Phys.  Rev.  Bl,  4005  (1970). 


A-4 


